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Abstract 

We present a density-functional theory (DFT) approach to the study of the 
phase diagram of the maximum density droplet (MDD) in two-dimensional 
quantum dots in a magnetic field. Within the lowest Landau level (LLL) ap- 
proximation, analytical expressions are derived for the values of the parame- 
ters N (number of electrons) and B (magnetic field) at which the transition 
from the MDD to a "reconstructed" phase takes place. The results are then 
compared with those of full Kohn-Sham calculations, giving thus information 
about both correlation and Landau level mixing effects. Our results are also 
contrasted with those of Hartree-Fock (HF) calculations, showing that DFT 
predicts a more compact reconstructed edge, which is closer to the result of 

exact diagonalizations in the LLL. 

PACS numbers: 73.23.Ps, 71.15.Mb, 73.40.Hm 
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Two-dimensional quantum dot systems, at high magnetic fields, have been recently stud- 
ied by various authors.&i The theoretical interest in these systems arises largely from the 
fact that they provide a few-electron realization of physical states which, in the macroscopic 
limit, are responsible for the occurrence of the quantum Hall effect. i@ The simplest example 
of such a state is the so-called maximum density droplet (MDD) — which, in the limit of 
high magnetic field, can be written as a Slater determinant of lowest Landau level orbitals 
with angular momenta 0, 1...N — 1, where N is the number of electrons.il In the limit of 
N — > oo this coincides with the incompressible state of the quantum Hall effect at filling 
factor v — 1. Because, within the LLL, the MDD is the only iV-electron state of angular 
momentum N(N — l)/2, (and there is none with lower angular momentum) it follows that 
it must be an exact eigenstate of the Hamiltonian 
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if the small Coulomb coupling between different Landau levels is neglected. Here Uq is 

the frequency of the external parabolic potential, Aj is the external vector potential, k the 

dielectric constant, m* the electron effective mass, /!# the Bohr magneton, g* the effective 

g-factor for the Zeeman splitting, and <jj the spin component along the axis perpendicular 

to the plane of the electrons. The question is whether this exact eigenstate (or rather 

its continuation to finite magnetic field) can actually be the ground-state of the quantum 

dot, in some range of magnetic fields. The basic physics is simple: if the magnetic field is 

too large, the MDD cannot be the ground state, because the compact arrangement of the 

electrons costs too much electrostatic energy: the electrostatic stress is released through a 

rearrangement of the electrons leading to a state of higher angular momentum. If, on the 

other hand, the magnetic field is too weak, the confinement energy will cause the external 

electrons in the MDD to be transferred to the center of the quantum dot, even though, in 

so doing, a higher Landau level becomes populated at the center of the dot. The conclusion 

of these arguments is that there will exist, at most, a "window" of magnetic fields in which 



the MDD is stable. The "window" shrinks with increasing electron number N and it closes 
up completely at a critical value N c of N of the order of 100. Note that this is not in 
contradiction with the existence of incompressible phases in the macroscopic limit: it is only 
telling us that such phases will have compressible edges. 

The problem of determining quantitatively the region of stability of the MDD has been 
studied both theoretically and experiment allyJl~i There exists a disagreement between the 
experimental results and theoretical predictions regarding the window of values of magnetic 
field for which the MDD is the ground state.! Correlation effects have been indicated as a 
possible cause for this disagreement since they were not accounted for in the theoretical anal- 
ysis, which was based on the Hartree-Fock approximation. The importance of correlations 
has been demonstrated for the case of small quantum dots.l 

Here we present an analytic treatment based on density functional theory, which includes 
both exchange and correlation effects. We shall calculate the values of the magnetic field 
at which the transition from the MDD to a new phase takes place, as well as the angular 
momentum of the new phase, supposed to lie entirely within the lowest Landau level. Within 
this treatment the transition from one state to the other can be entirely described by means 
of a single dimensionless parameter a = m*uj 2 k\ 3 /e 2 giving the strength of the parabolic 
potential in terms of ratio between the confining energy and the electrostatic energy existing 



at the typical length scale of a magnetic length A = J he/ (eB). We also perform a numerical 
evaluation, based on the solution of the Kohn-Sham equation, where all Landau levels are 
included, showing that their inclusion shrinks the magnetic field window of stability of the 
MDD. As a byproduct of our approach we determine the maximum number of electrons for 
which the MDD can be the ground state. 

Density functional theory has already been applied successfully to systems in the presence 
of a magnetic fieldJ^HU thus establishing it as a useful tool for studying such systems. The 
total energy of a quantum dot is 
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In Eq. (0) we have assumed the system is within the lowest Landau level, thus omitting the 
constant kinetic energy term tkv c /2, where u c = eB/(m*c) is the cyclotron frequency. V(r) 
represents the parabolic confining potential 

V(v) = \m^y = |r 2 -^, (3) 

The dimensionless parameter a gives the strength of the parabolic potential. ^ c [n] is the 
exchange-correlation energy functional. Let now tim dd(y) be the density of the MDD, and 
n r (r) the density of the reconstructed edge obtained immediately after the transition from 
the MDD takes place for example by increase of the magnetic field. In the lowest Landau 
level picture the MDD is obtained by filling the orbitals with angular momentum from I = 0, 
to I = N — 1, the system being fully spin polarized. The reconstructed edge can be viewed 
as being generated from the the MDD by removing one electron from an orbital <ph with 
angular momentum < h < iV — 1, and putting it in the single-particle orbital 0jv with 
angular momentum N. Therefore 

n r (r) = n MDD (r) - (^(r)] 2 + |<Mr)| 2 . (4) 

The transition from one state to the other for a given number of electrons is totally deter- 
mined by the value of the parameter a. Its critical value a% is obtained by solving 

E[n MDD , cti, N) = E[n r , a u N]. (5) 

Using Eq. (0) and @, the energy difference between the two states can be written as 

E[n MDD , a, N] - E[n r , a, N] = 
dv V{v) [|<Mr)| 2 - |<Mr)| 2 l + AE H + AE XC . (6) 



We are now going to evaluate separately each term on the right hand side of Eq. (Q). The 
first term is easily obtained from the second moment of the Landau orbitals 

Rf = |drr 2 |<Mr)| 2 = 2A 2 (|Z| + l). (7) 



The change in Hartree energy AEh is given by 
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In order to proceed we approximate the MDD as a system with uniform density, having 



the N electrons in a disk of radius Rmdd = 
by a disk of radius R of uniform charge uq is 

r m > e ( R - r ') 



f 2N\. The electrostatic potential generated 
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where and E(x) are the complete elliptic integrals of the first and second kind respec- 

tively. The densities associated with the orbitals </»^(r) and </>jv( r ) are treated as properly 
normalized delta functions |^(r)| 2 ~ 5(r — Ri)/(2irRi), in the terms of Eq. (|8]) which in- 
volve the interaction of any of these two orbitals with the MDD and with each other. The 
self- interaction of the orbitals 0;'s is treated as that of rings whose electrostatic energy is 
e 2 /{2kCi) with the capacitance of the ring C\ = tt Ri/ [In (2Ri/X)]. Using these approxima- 
tions in (If) we obtain 
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The variation of exchange-correlation energy between the two states can be evaluated within 
a local density approximation 

AE XC = jdr e xc (l) 9(R MDD - r) (n MDD (r) - n r (r)) 

2 

= jdr e xc (l) 9{R M dd - r) (|0,(r)| 2 - |<Mr)| 2 ) » e xc (l) = 6^(1) (11) 

where we have neglected variations in the tail of the MDD and we have again replaced the 
orbital densities with S— functions. e xc (l) is the exchange-correlation energy per particle of 
a uniform electron gas for filling factor v = 2im{r)l 2 = 1. From Eqs. (0), (|10|), and ([□]), we 
arrive at 
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By equating the left-hand side of Eq. ( |12"1) to zero, namely by looking at the transition from 
the MDD to the reconstructed edge, one gets the value of a = a±(N; h) which characterizes 
this transition. The reconstructed phase occurs when a < ai(N; h), for a given value of h. 
Therefore the position of the first hole in the reconstructed edge is obtained by maximizing 
ai(N; h) with respect to h. This permits also to derive the maximum possible value c**(iV) 
for which a transition from the MDD to the reconstructed edge is possible. If N is sufficiently 
large and h ~ N, we have 
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Differentiation with respect to h gives 
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with 2-^1(0, c; z) the hypergeometric function. Under the hypothesis — which we shall 
show is valid — of h ~ N, the ratio (h + 1)/N is a number smaller than one, but close 
to unity Therefore — in order to find an expression for the change of angular momentum 
associated with the reconstruction of the edge — it makes sense to use an expansion of the 
hypergeometric function in terms of 1 — z, and retain the lowest order terms. FromB 

E[x] wl + (l-i)[oi- 61 In (1 - ar)]; a x = 0.44325141463, b x = 0.24998368310, (15) 

and recalling the representationlll 
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for the hypergeometric function, valid whenever c = a + b + /, I — 0, 1, 2, . . . , we obtain 
retaining only the lowest order term in (1 — z), 
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In Eq. 



k n = ^(n + 1) + ^(n + 1 + 1) - ip(a + n + I) - tp(b + n + I). 



(18) 



ip is the logarithmic derivative of the gamma function, while (a) n = T(a + n)/T(a), is 
Pochhammer's symbol. Application of the condition dai/dh = 0, to Eq. ( |i~7|) results in 
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where we have used 4 bi ~ 1. The value of e sc (l) = — 0.7015e 2 / (fcZ), gives N — h w 3.115 V^V- 
The same scaling behavior was found numerically by Oaknin et a/.Jilwith the only difference 
in the value of the prefactor, being 2 in their case. The difference between our and their 
results might be a consequence of the fact that the latter have been obtained for systems 
smaller than those in which the asymptotic behavior of Eq. fll9|) applies. We also notice 
that introducing the result of Eq. (^) in Eq. (0), the latter agrees with the behavior for 
a\ predicted by Chamon and Wen,! except for the value of the numerical factors. The 
difference in the prefactors is due to the fact that while in our case we obtain it from an 
asymptotic behavior, Chamon and Wen determined it in terms of a fitting procedure. In 
Fig. [1] we present the increase of angular momentum with respect to the MDD after the 
edge reconstruction takes place, scaled by the square root of the number of electrons, as a 
function of the number of electrons in the dot, obtained by numerically finding the maximum 
with respect to h of cti(iV; h) as obtained by equating to zero the left-hand side of Eq. (|i2|). 
The asymptotic value of Eq. ([H]) for large dots — represented by a dashed line in Fig. [l 
- is approached very slowly, basically for values much larger than those considered in the 
figure. 

The other boundary of the MDD region of the phase diagram is derived in an analogous 
fashion. Here the new phase is described by flipping the spin of one electron in the MDD 
and putting it into a state with angular momentum < I < N — 2. The transition between 
the two phases occurs when 



In Eq. (|20|) , the last term accounts for the change in exchange-correlation energy, while 




(20) 
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the remaining ones are purely electrostatic, and represent the energy required for moving a 
charge distribution located around Rn-i to a neighborhood of R[. We see that this region 
of the phase diagram is dominated by electrostatics, exchange-correlation effect being of 
higher order in l/y/N. The expression for Eq. ( p0|) has to be minimized with respect to /, 
in order to obtain the first configuration that gives the ground state when the MDD is lost 
for increase of a. Considering the terms with leading order in N, we arrive at 
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Eq. (|2l|) is monotonically increasing with I. It follows then that the angular momentum that 
minimizes Eq. ([H]) is I — 0. Therefore 

1 ( 2\/2N 1 
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(22) 



From the previous analysis we conclude that the MDD is a ground state whenever a\ < 
a < ol\. Since the expressions for a\ and a\ depend only on the number of electrons N, the 
maximum value N c for which the MDD is a ground state for some values of a is obtained 
from «2(A^ C ) = al(N c ). By solving the previous equation by means of (0) and (p2|), we get 
N c = 222. If instead we use the expressions for the a*'s as obtained by equating to zero (15) 



and from (|20f) , finding then their maximum with respect to h, and minimum with respect 
to I respectively, we obtain N c = 160. 

The above discussion was limited to the lowest Landau level. We now turn to considering 
the effects coming from the inclusion of higher Landau levels. In this case the full Kohn-Sham 
equations for the quantum dot must be solved.0 The results for the B-N phase diagram of 
the MDD are presented in Fig. [| The most prominent feature (in comparison to the LLL 
approximation) is a narrowing of the window of magnetic fields for which the MDD is the 
ground state. In particular, N c ~ 28, when ujq = 4meV. 
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Therefore we see that Landau level mixing is essential to an accurate determination of 
the correct values at which the transitions take place. Correlation effects are manifested 
by giving a more compact quantum dot with respect to what is given by Hartree-Fock 
theory, which includes only exchange. Moreover, the values of the angular momentum of 
the reconstructed edge as predicted by DFT are in better agreement than HF with exact 
diagonalization calculations. This is represented in Fig. ^ where we give the values of 
the energy of a MDD with a hole in it (i.e., a MDD on the verge of reconstruction), as 
a function of the position of the hole, for three different values of the magnetic fields, as 
evaluated by numerical diagonalization within the LLL, by DFT, and by HF. AM is the 
increase of angular momentum for "one-hole" states with respect to the angular momentum 
of the MDD. The fact that DFT departs from the exact diagonalization results for large 
AM is not surprising, since those states correspond to excited states, where DFT is not 
applicable. Of course these states do not affect the B-N phase diagram for the MDD. 

In conclusion, we have shown that correlation effects which are incorporated in DFT 
give rise to a more compact reconstructed edge than the one predicted by HF theory. These 
results are in better agreement with those of exact diagonalization studies. A simplified 
model reproduces the main physical traits involved in the phase diagram of the MDD, 
which is determined, within the model, by means of analytical expressions. 
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V. acknowledges financial support from NSF under grant No. DMR-9403908. 
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FIG. 1. Change of angular momentum N — h in going from the MDD to the reconstructed edge 
scaled by y/N, as obtained by maximizing a\(N;h) with respect to h. The dashed line represents 
the asymptotic value of Eq. 
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FIG. 2. B-N phase diagram of the MDD for ujq = 4eV obtained solving the Kohn-Sham equa- 
tions as in Ref. [7^. 
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FIG. 3. Evolution of the energy for the MDD and the "one-hole" states as a function of the 
magnetic field for a 20 electron quantum dot. ujq = 3meV. 
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